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We consider the effects of temporal delay in a neural feedback system with excitation and inhibition. The 
topology of our model system reflects the anatomy of the avian isthmic circuitry, a feedback structure found in all 
classes of vertebrates. We show that the system is capable of performing a 'winner-take-all' selection rule for 
certain combinations of excitatory and inhibitory feedback. In particular, we show that when the time delays 
are sufficiently large a system with local inhibition and global excitation can function as a 'winner-take-all' 
network and exhibit oscillatory dynamics. We demonstrate how the origin of the oscillations can be attributed 
to the finite delays through a linear stability analysis. 



I. INTRODUCTION 



In order to identify and react to behaviorally relevant ob- 
jects in their visual environment, animals must be able to 
rapidly locate the positions of these objects in visual space. 
This ability to select and orient towards the most salient 
part in a visual scene that may be cluttered with other, 
for the animal's survival less relevant objects, has evolu- 
tionary significance, as it permits the organ ism to detect 
quickly possible prey, predators, and mates ('i tti and Kochi 
I2OOI). In standard models of selective visual attention, the 
stimulus is encoded in a 'saliency map' that topographi- 
cally represents the conspicuity of the stimulus over the vi- 
sual scene. The most salient location is then chosen by a 
'winner-take-air (WTA) n etwork, i.e., by a neurall y imple- 
mented maximum detector jKoch and UIIman[|l985h . In neu- 
ronal network models, these WTA n etworks are often r e alized 
as networks w ith lateral inhibition ( Amari and Arbibl ll977t 
ICouItrip et all [1992: Kaski an d Ko honen. 1994 ,). global inhi- 
bition ( Indiveri and Delbriickl 2002h or l ocal ex citation and 
long distance inhibition ( Standage et al.l l2005b . After the 
most active location, i.e., the 'winner', in the saliency map 
has been chosen, attention should not, however, continue to 
be focused onto it. One way of allowing attention to shift, is 
to transiently inhibit neurons in the saliency map that corre- 
spond to the currentl y attended locat ion, a strategy known as 
'inhibition of return' (lltti et aI.L[l998l) . 



The homolog of the mammalian superior coUiculus in 
non-mammaUan vertebrates is the optic tectum (TeO). It 
is critically involved in localizing visual objects and in 
the preparation of orienting responses to wards these objects 
(ISparksl Il986l; IWurtz and Albanol Il980l) . In all classes of 
vertebrates, the TeO is reciprocally connected with the nu- 
cleus isthmi (N I), which is hom ologous to the parabigeminal 
nucleus (Diam ond et all [T992I) in mammals. In the avian vi- 
sual pathway, the NI consists of three subnuclei: the nucleus 
pars parvocellularis (Ipc), the nucleus pars ma gnocellularis 
(Imc) , and t he nucleus pars semilunaris (SLu) (IWang et all 
I2OO6 , I2OO4I) . In both Ipc and Imc the projection from the 
tectum is topographically organized such that the retinotopic 
map is preserved in both nuclei, with the projectio n to the 
Imc being somewhat coarser than for the Ipc (Wan g et"all 
I2OO4I) . In contrast, the isthmic projections back to the TeO 
are very different for Ipc and Imc. Ipc neurons project back 



to the TeO in a highly precise homotopic manner, i.e., the 
axons of each Ipc neuron terminate in that part of the op- 
tic tec tum from which their visual inputs come (Wang etall 
l2006h . Imc, on the other hand, has two populations of neu- 
rons, which both make h eterotopic projectio ns but only to the 
TeO or Ipc, respectively ('Wang et"alll2004l ). The three-nuclei 
circuitry consisting of TeO, Ipc, and Imc is shown in Fig. [T] 
(iMarm et all l2005b . Due to latencies arising from synaptic 
process and the spatial separation of the nuclei, the coupling 
between TeO and NI cannot be considered instantaneous . 
Rather, finite temporal delays exist (IWang etalll2006l |2004 . 
Furthermore, delays can arise from the dy namical proper- 
ties of the systems involved. For instance, lAndersen et al.l 
(1964) report stimulus-dependent onset latency of recurrent 
inhibition in the c at hippocampus, and these findi ngs were 
later explained by iHauptmann and MackevI (l2003h . It has 
been known for some time that temporal delays can cause an 
otherwise stable system to oscillat e (Coleman and Renn ingej 
119761: iHadeler and TomiukL Il977t Ian der Heidenl [19791) and 
may lead to bifurcatio n scenarios re sulting in chaotic dynam- 
ics dSchanz and Pelst er, 2003; Wis chert et all Il994 . There- 
fore, finite delays are an essential pr operty of any realistic 
model of a neuron population ( Miltoniri996h . 

The synaptic effect of the recurrent projections from the 
Ipc and Imc cells onto their target cells is less well under- 
stood than their anatomical organization. The available evi- 
dence suggests that Ipc neurons are cholinergic, whereas Imc 
neurons have been shown to express gamma-aminobutyric 
acid (GABA) as their main neurotransmitter Thus, accord- 
ing to the usual role of acetylcholine and GABA, one might 
speculate that Ipc and Imc neurons mediate ex citation and 
inhibit ion onto their target cells, respectively. iMarin etal] 
(l2005b posit that "the three-nuclei circuitry [...] may con- 
stitute a winner-take-all network iKoch and UIImanLTT985l) 
in which local visual inputs to the Ipc are augmented by 
the re-entrant loop among tectal and Ipc neurons, combined 
with broad inhibition of the rest of the Ipc by Imc neu- 
rons." This argument seems imme diately plaus i ble, h o wever , 
electrophysiological e xperiments ("Wan g et all 1 19951 l2000t 
IWang and Frost , 1 1991,) suggest that the synaptic effects of 
Ipc and Imc are actually converse to this scenario and that 
the Ipc mediates inhibition whereas the Imc has an excita- 
tory effect. Given the anatomical organization of the re- 
current projections from Ipc and Imc, it is not fully in- 
tuitive how the system could function as a WTA network 
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FIG. 2: Diagrammatic representation of our model for the istlimo- 
tectal feedback loop. Neurons in the TeO and Ipc, which project 
topographically, are modeled as individual units. Due to their diffuse 
projections, Imc neurons are combined to form a feedback kernel 
(cf. Tab. IB. 



FIG. 1: Diagram of the isthmotectal feedback loop. The Ipc is recip- 
rocally connected with the TeO in a precise homotopic manner. Tec- 
tal neurons project topographically to the Ipc, and Ipc neurons project 
back to the corresponding tectal loci. Imc receives a coarser topo- 
graphic projection and projects back to the TeO and Ipc via widely 
ramifying terminal fields. Black represent s visually ac t ivated neural 
elements. Reprinted with permission from lMarin et al.l bOOSi) . 



when Imc i s excitatory and Ipc inhibitory. Nevertheless, 
IWangI (l2003b considers this possibility: "The positive and neg- 
ative feedback loops formed between the tectum and NI may 
work together in a winner-take-all network, so that the posi- 
tive feedback loop could provide a powerful augmentation of 
activated loci, while the negative feedback loop may strongly 
suppress the others [...]. For example, Imc could enhance the 
visual responses of tectal cells to target locations or stimulus 
features, while Ipc may suppress those to other locations or 
features in the visual field." The aim of this work is to investi- 
gate possible mechanisms for WTA selection in the isthmotec- 
tal feedback loop through a computational model. In this con- 



TABLE I: Components of the isthmotectal feedback loops and ab- 
breviations. We use the greek indices a, (3, and 7 to denote TeO, Ipc, 
and Imc, respectively. Furthermore, neurons are numbered such that 
indices 1 through A'^ refer to the TeO, A'^ + 1 through 2N refer to the 
Ipc, and the index 2N + 1 refers to the Imc. 



optic tectum nucleus isthmi pars nucleus isthmi pars 
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a 


P 
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1...N 


TV + 1 ... 27V 


27V + 1 



text, we do not refer to the term WTA in its most strict sense, 
which would imply that only the neuron with the strongest in- 
put exhibits a nonzero firing rate; rather, we speak of WTA 
behavior when the firing rates of neurons with weaker inputs 
are suppressed relative to those with stronger input. 

In Sect, ini we introduce our model of the isthmic system, 
and we analyze its response dynamics for different temporal 
delays and different combinations of excitation and inhibition 
in Sect.ini] In Sect.|IV] we compare the efficiency of WTA se- 
lection for these combinations. In Sect.|V] we employ a linear 
stability analysis to show how the oscillatory dynamics that 
arise in the system can be attributed to the increasing delays. 
In Sect. IVll we summarize our results. 



II. MODEL 

To explore the conjecture that the isthmotectal feedback 
loop functions as a WTA network, we consider a model 
system of coup l ed Hopfield neurons with tem poral delays 
(iHopfieldi IT981 iMarcus and Westervelj. 1 1989b . as described 
by ErmentroutI ( 19981) . In this model, the temporal evolution 
of the membrane potential of the ith neuron (taken from rest 
potential), Vi{t), is given by the first-order delay differential 
equation (DDE) 



dt 



V^{t) + yw^,rj{t ~ nj) + h{t) . ( 1 ) 



3 



Here, the membrane time constant for the ith neuron is de- 
noted by T^"^\ the synaptic connection weights for the pro- 
jection from the jth to the ith neuron are Wij, the temporal 
delay for this projection is r^j, is the firing rate for the jxh 
neuron and is linked to its voltage according to a nonlinear 
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firing rate function, 



(2) 



and Ii{t) denotes an external input to the ith neuron. To 
model the isthmic system, we assume that N tectal neurons 
are reciprocally coupled to N Ipc neurons and that the only 
neurons that receive external input are those in the TeO. Fur- 
thermore, due to the broad and heterotopic nature of the pro- 
jections from Imc, we combine the Imc neurons to a feed- 
back kernel, which then projects globally to both TeO and 
Ipc. The topological structure of our model is depicted in 
Fig.|2] To simplify our model, we make the following assump- 
tions: The synaptic weights for the projections TeO— >Ipc, 
TeO^Imc, Ipc— > TeO, Imc^TeO, Imc^Ipc, are the same 
for all neurons in each of these groups, and we denote them 
by wpa, w.ya, Waj3, Wa-y, Wfj^, respectively; all membrane 
time constants are identical, r^^™^ — r'"^^ for all i, and we 
rescale time such that r^'"^ = 1; all delays are identical, 
Tij = r; all firing rate functions are identical ^^(V,) = 
S{Vj). Furthermore, we number our neurons such that the 
indices i = 1, 2, . . . , iV refer to tectal neurons, the indices 
i = A^ + l, iV + 2, 2N refer to Ipc neurons, and the in- 
dex i = 2N + 1 refers to the Imc kernel (cf. Tab.|l|. Then, the 
dynamics of our system are described by the 2N + 1 DDEs: 



dt 

+Wa^r2N+i{t -t) + Ii{t), i = 1,2, ... ,7V 
dV,{t) 
dt 

+wp^r2N+iit -t), i^N+l,N + 2,...,2N. 

dV2N+l{t) ^ 



dt 



= -Vi{t) + Wapr,+N{t - t) 

{t-T)+h{t), i=l, 
Vi{t) + WpaTi-Nit - t) 

(t-r), ^ = iv + l,iv^ 

N 

a. 

Y.n{t-T). (3) 



i=i 



For the firing rate function we choose the piecewise linear 
function. 







for V, < Vt, 



S{V,) = { a{y, - Vt) for Vt < V, <Vt + S^^^/a, 

max 

for S'max/ a + VT<Vj. 

(4) 

Finally, we make the simplifying assumptions \wpa\ = 
N Iwjal = \wai3\ = Iwajl ~ Iwp-y] = l/fl and Vt — 0. The 
signs of the synaptic weights determine whether a projection 
is excitatory or inhibitory and since the tectal cells me diate 
excita t ion, w e have wpa, w-ia > 0. iMarm etaL I (120051) and 
IWand (l2003h both discuss scenarios in which WTA behavior 
arises from an interplay of excitation and inhibition in the isth- 
motectal feedback loop, and we are therefore interested in the 
cases where Ipc and Imc have adversary effects onto the TeO. 
Thus, four cases remain to be studied, which can be character- 
ized according to the signs of {wap, Wa-y, wp^) as (+, — , — ), 
(+, — , +), (— , +, +), and (— , +, — ). In our model, the first 
two of these cases correspond to global inhibition and local 
excitation of tectal cells through feedback, whereas the latter 
two correspond to global excitation and local inhibition of the 
cells in the TeO. 



ni. RESPONSE DYNAMICS 

In the following, we investigate the dynamical response 
behavior of a network consisting of = 200 (initially 
quiescent) neurons to a static stimulus. We choose an in- 
put consisting of five superimposed Gaussians with peaks at 
i = 20, 60, 100, 140, and 180, and peak values of 0.75, 0.5, 
0.45, 0.4, and 0.35, respectively. The (normalized) stimulus is 
shown, e.g., in Fig.[3](a). In particular, we are interested in the 
firing rates of those neurons whose positions coincide with the 
peaks in the stimulus, and in order to abbreviate our notation 
we denote their indices as a, h, c, d, and e according to de- 
scending strength of their respective inputs. 



A. Global inhibition, local excitation 

We first consider the cases (+,—,—) and (+,—,+). In 
this situation, our network is similar to the circuit considered 
in (Indiveri and Delbriick, 2002). Therefore, we expect that it 
can perform a reasonably accurate WTA selection. Figure [3] 
shows the firing rate dynamics in response to the static input. 
The undelayed case is shown in Figs. [3] (a) and (b), whereas 
Figs.[3](c) and (d) show the dynamics that result when r = 2. 
To compare differences between firing rates of tectal neurons 
we consider the contrast measure 



Ti + r,- 



(5) 



In particular, we are interested in the contrasts Cab, Cac, Cad, 
and Cae, whose temporal dynamics are depicted in the second 
plots of Figs. [3](b) and (d). From Fig. [3] we see that in the 
case (+,—,—) the weakest inputs are suppressed efficiently, 
while the neurons receiving the strongest input are driven to- 
wards maximum firing. Inputs of intermediate strength (e.g., 
the one received by neuron b) are not suppressed. Thus, in the 
configuration (+,—,—) the system can perform a WTA se- 
lection, but not with very good 'resolution' . Furthermore, by 
comparing Figs.[3](a) and (b) with Figs.[3](c) and (d), we see 
that the temporal delay in the system has only little effect on 
its efficiency as a WTA selector. The main effect of the delay 
is that it causes the system to evolve on a longer time scale, 
i.e., the steady state is not reached as fast as in the undelayed 
case. 

For the case (+,—,+) it turns out that the inhibition in the 
system is insufficient to compensate for the positive feedback 
in the recurrent coupling between TeO and Ipc, and even neu- 
rons that receive only weak inputs are driven towards maxi- 
mum firing. Thus, in the configuration (+, — , +) our model 
system does not function as a WTA network. 



B. Global excitation, local inhibition 

Now we consider the cases (—,+,+) and (—,+,—), which 
correspond to global excitation of the TeO through the Imc 
and local inhibition from the Ipc. They thus correspond to 
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FIG. 3: WTA behavior and network dynamics for the case (+, — , — ), i.e., global inhibition and local excitation, (a) and (b) are for the 
undelayed case, r = 0, while in (c) and (d) the delay is r = 2. The dots in the snapshots in (a) and (c) show the normalized response of tectal 
cells, i.e., the quantity fi = ri{t)/rina.x{t), where rmax(t) denotes the maximum of all tectal firing rates at time t, while the solid line shows 
the normalized input, i.e., the quantity li/Imax- In both (b) and (d), the first plot shows the firing rate dynamics of the neurons a through e in 
the TeO, the second plot depicts the temporal evolution of the response contrast of tectal neurons b through e when compared with neuron a, 
the third plot shows the activity of the Ipc neurons receiving input from the tectal neurons a through e, and the fourth plot shows the firing rate 
of the Imc. 



the scenario described bv lWand (l2003h and it does not seem 
immediately intuitive how WTA behavior could result in this 
configuration. Indeed, for the case (— , +, +), when there are 
no delays no WTA selection occurs. The response dynam- 
ics for this case are shown in Figs. |4] (a) and (b). The con- 
trast in the firing rate response of neurons receiving inputs 
of different strengths is nearly identical to the contrast of the 
respective inputs during all phases of the system's temporal 
evolution. Thus, neither are weaker inputs suppressed nor are 
stronger inputs augmented. When we introduce delay into the 
system, however, its behavior changes drastically. Figs.|4](c) 
and (d) show the response dynamics for the case r = 2. Note 
that the neurons' firing rates, as well as the contrasts between 
responses now exhibit oscillatory behavior The system can 
perform a WTA selection with reasonable accuracy, but only 
transiently, i.e., only during certain phases of its temporal evo- 
lution. As a matter of fact, a phase of best WTA selectivity is 
preceded and followed by phases where the response contrast 
is even lower than that of the input. 

Since standard models of selective attention usually require 
that the most salient stimulus not be a permanent 'winner', but 
rather that it be suppressed once attention has been directed 
to it, the dynamical evolution of the response is an important 
characteristic. Consequently, it may actually be a desirable 



feature of a WTA network to only determine the 'winner' tran- 
siently. 

The case (— , +, — ) leads to similar results as in the case 
(+, — , +). When Ipc neurons are inhibited by the Imc, they 
cannot provide sufficient negative feedback to the TeO in or- 
der to prevent tectal neurons from being saturated through the 
positive feedback between TeO and Imc. Thus, in the con- 
figuration (— , +, — ) our model system does not function as a 
WTA network. 



IV. COMPARISON OF WTA SELECTIVITY 

In order to quantify the performance of our model system 
as a WTA network, we consider the maximum in the response 
contrast between neurons i and j, normalized to the contrast 
between the (constant) input / to neurons i and j during the 
first 30 membrane time constants of the system's temporal 
evolution: 



I, 



\Ii — Ij \ 0<i<30 



(6) 



This quantity is shown for the pairs (a, b), (a, c), (a, d), and 
(a, e) for the cases (+, — , — ) and (— , +, +) and for differ- 
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FIG. 4: WTA behavior and network dynamics for the case ( — , +, +), i.e., global excitation and local inhibition. Same as Fig. [3] but for 
inversed signs of u)q,/3, Wa-,, and wfj-,. 



ent values of the time delay in Fig.|5] As we expect from the 
results presented in Sec. [Till in the case of global inhibition 
and local excitation, (+, — , — ), the system's performance as 
a WTA network, measured by the value of depends only 
little on the time delay. Furthermore, we see that the system 
is efficient in suppressing weak inputs, whereas the response 
contrast for inputs of intermediate strength is less enhanced. 
The performance for the case (— , +, +), on the other hand, 
depends strongly on the temporal delay. Furthermore, the ra- 
tio between the maximal response contrast and the input con- 
trast is comparable for weak and intermediate inputs. When 
the delay is sufficiently large, the model system thus exhibits 
a better WTA 'resolution' in the case of global excitation and 
local inhibition than for the inverse scenario, albeit only tran- 
siently. 

We have also investigated the role of parametric disorder in 
the system and find that it does not change our results qualita- 
tively. For instance, when the projection latencies are drawn 
randomly from a normal distribution with mean at r = 2 
and standard deviation of 0.2 and the synaptic weights are 
disorderd using normal distributions with means at the de- 
fault values and standard deviations of 10% of these val- 
ues, the maximum contrast as measured by is compa- 
rable to the case without disorder. Simulating ten different 
samplings of randomized delays and synaptic weights for the 
case (-,+,+), we obtain the result {Cab, Cac, Cad, Cae) = 
(2.28 ± 0.25, 2.50 ± 0.30, 2.06 ± 0.13, 2.34 ± 0.12) (re- 
sults are mean ± standarad error of the mean), which is 



to be compared with the values for the undisordered case 

{Cab, Cac, Cad, Cae) = (1-83, 1.89, 2.00, 2.06). 



V. LINEAR STABILITY ANALYSIS 

We now aim to understand the origin of the delay-induced 
oscillatory dynamics in the case of global excitation and local 
inhibition through a stability analysis of the model system. To 
this end, we make the following ansatz, which, a posteriori, 
turns out to be correct. We assume that for the chosen input 
the system of DDEs ([3]) possesses a stationary point Vi [t) = 
Vi with Q <Vi < S'max/a for 1 < j < 2iV + 1, and we can 
thus replace the voltages in the system (O according to the 
linear part of the firing rate function S. The stationary point 
is then obtained by solving the equation 



V = aW\ + I . 



(7) 



Here, V and I are 2N + 1 column vectors (only the first N 
entries of I are nonvanishing) and W is a (2iV+ 1) x {2N + 1) 
matrix of the form 



W = 



OnxN 
Wja 1 1 X TV 



OlxN 



(8) 



We find that for the case (— , +, +), the matrix 1 — aW is 
invertible and that the solution for the stationary point V = 
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FIG. 5: Performance of the model system as a WTA network. 
The maximal response contrast normalized to the input contrast is 
shown as a function of the delay r. Open symbols are for the case 
(+, — , — ), and filled symbols are for the case ( — , +, +). The con- 
trast Cij is shown for i = a and j — b (circles), j — c (triangles), 
j — d (squares), and j — e (diamonds). 



(1 — aW)~^I does indeed permit us to linearize the system 
(O. In the case (+,-,-), however, it turns out that this solution 
yields values that lie outside of the linear regime of the firing 
rate function, and the linearization Q is therefore not valid. 
Next, we analyze the stability of the stationary point for the 
case (— , +, +) by making the ansatz V(t) = V + ce"^*, which 
leads to the equation [ae^^'^W — (1 + A)!] c = 0. In order 
to determine the conditions for a nontrivial solution to this 
equation to exist, we must solve the characteristic equation 
for the matrix M — ae~^'^W — 1, i.e., we have to solve 



det(M- Al) 



M 
u 



V 

-(1 + A) 



(9) 



where we have introduced the abbreviations 



M 



I -{1 + X)InxN aWafje ^''i-NxN 



^'^^Nxl 



(10) 

(11) 

(12) 



Solving (|9]l is facilitated by applying the identity d Pasolovl 
11994 det(M) = - det(M) [l + A - uM~iv] . The inverse 
of M is given by 



= -- 



1 



-2\t 



{1 + \)Inxn awape ^''Inxn 



(1 + A)2 - a^wapwpae 

)11 

\aWi3ae^^'^tNxN (l + A)lArxJV 

and we thus obtain 



(13) 



FIG. 6: Eigenvalues of M as a function of the delay r. Real parts 
of eigenvalues are shown by solid lines, imaginary parts by dashed 
lines. The upper panel shows real and imaginary part of A3, and the 
lower panel shows real and imaginary of Ai- and A2-. 



det(M- Al 



-At 



,(1 + A)] 



[aWapWp-fC 

X Na^w^a^^'^^'^ [(1 + A)^ — a^WaplVpaS-'^ 
- (1 + A) [(1 + \f - a'w^pwp^e-^^^] 



2\t 



N-1 



(14) 



For the c ase (—,+,+) the characteristic equation simplifies to 



[1 + e3^^(l + A)3] + (1 + A)2] = 0. Its solutions, 

the eigenvalues of M, are given by 



Ai± = -1 + -VF(±ire^) , 

T 



A2± = -1 + -W 

T 



A3 = -1 + -W{-eW). 

T 




(15) 
(16) 
(17) 



Here, W{z) is the inverse function to W^^(z) — ze^, which 
is usually called the Lambert ly-Function. Figure|6]shows the 
real and imaginary parts of the eigenvalues of M as a function 
of the delay r. For all values of r there is no eigenvalue with 
a positive real part and the fixed point V is thus stable for 
arbitrary delays. However, the real parts of the eigenvalues 
tend to zero faster than their respective imaginary parts as the 
delay is increased. Therefore, with increasing delay, the relax- 
ation time for the system's return to the stationary point grows 
more rapidly than the time scale for oscillations. Thus, with 
increasing delay, the system will spkal toward the fixed point, 
explaining the observed oscillatory behavior. 



VI. SUMMARY AND DISCUSSION 

We have investigated the circumstances under which the 
isthmic system can function as a WTA network. We have con- 
structed a rate-model of the isthmotectal feedback loop and 
have analyzed the temporal evolution of the model system in 
response to a static stimulus. We have shown that time de- 
lays can be crucial to the dynamical behavior of the system. 



7 



In particular, delay-induced oscillations can lead to transient 
WTA selection in our model. Finally, we have performed a 
linear stability analysis explaining the origin of the oscillatory 
behavior 

It has been conjectured for a long time that the isthmotectal 
feedback loop constitutes a WTA network. Our results show 
that the isthmic circuitry is indeed set up to perform such a 
selection rule. In the case where global inhibition and local 
excitation are present in the system, this result is quite intu- 
itive. However, a network with global excitation and local 
inhibition might not appear to be well-suited as WTA selec- 
tor. Yet, precisely such a scenario was discussed in the lit- 
erature. Temporal delays can be crucial for the behavior of a 
dynamical system, and, as we have shown in our investigation, 
they are particularly important for the case of global excitation 
and local inhibition, as they induce transient WTA behavior in 
the network. Transmission and synaptic delays for the projec- 
tions between Ipc and TeO are estimated to be around 15ms 
(fNetzel et al...,2006.) . whereas membrane time constants in the 



Ipc may be as short as a few milliseconds (J. Shao, personal 
communication), which is within t he ran ge of typical neuronal 
membrane time constants (lKochLll999h . When the synaptic 
and transmission delays are of the same order of magnitude 
as the membrane time constants involved, the degree to which 
our model for the isthmic system functions as a WTA network, 
can depend crucially on the delays (cf. Fig. |5]l. Therefore, 
temporal delays should not be neglected when the neuronal 
dynamics of the isthmotectal feedback loop are assessed and 
its potential for WTA selection is discussed. 
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